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The problem of non-perturbative description of stationary flames with arbitrary gas expansion 
■ is considered. On the basis of the Thomson circulation theorem an implicit integral of the flow 

equations is constructed. With the help of this integral, a simple explicit expression for the vortex 
mode of the burnt gas flow near the flame front is obtained. Furthermore, a dispersion relation 
£SJ ■ for the potential mode at the flame front is written down, thus reducing the initial system of bulk 

equations and jump conditions for the flow variables to a set of integro-differential equations for the 
\ flame front position and the flow velocity at the front. The developed approach is applied to the 

Qcase of thin flames. Finally, an asymptotic expansion of the derived equations is carried out in the 
case 8 — > 1 where 8 is the gas expansion coefficient, and a single equation for the front position is 
£ — ' obtained in the second post-Sivashinsky approximation. It is demonstrated, in particular, how the 

£N| ■ well-known problem of correct normalization of the front velocity is resolved in the new approach. 

It is verified also that in the first post-Sivashinsky approximation, the new equation reduces to the 
Sivashinsky-Clavin equation corrected according to Cambray and Joulin. Analytical solutions of 
the derived equations are found, and compared with the results of numerical simulations. 
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.£h '. I. INTRODUCTION 

>y 

The process of flame propagation presents an extremely complicated mathematical problem. The governing equa- 
O |' tions include the nonlinear flow equations for the fuel and the products of combustion, as well as the transfer equations 
l— ~~ ', governing the heat conduction and species diffusion inside the flame front. Fortunately, in practice, an inner flame 
' scale determined by the latter processes is large compared to the flame front thickness, implying that the flame can be 
considered as a gasdynamic discontinuity. The initial problem is reduced thereby to a purely hydrodynamic problem of 
t^J- ■ determining the propagation of a surface of discontinuity in an incompressible fluid, the laws of this propagation given 
by the usual Navier-Stokes and continuity equations complemented by the jump conditions at the surface, expressing 
the mass and momentum conservation across the flame front. The asymptotic methods 1-3 allow one to express these 
conditions in the form of a power series with respect to the small flame front thickness. 

Despite this considerable progress, however, a closed theoretical description of the flame propagation is still lacking. 
What is meant by the term "closed description" here is the description of flame dynamics as dynamics of the flow 
variables on the flame front surface. Reduction of the system of bulk equations and jump conditions, mentioned above, 
to this "surface dynamics" implies solving the flow equations for the fuel and the combustion products, satisfying 
given boundary conditions and the jump conditions at the flame front, and has only been carried out asymptotically 
for the case 6 — ► 1, where 9 is the gas expansion coefficient defined as the ratio of the fuel density and the density of 
burnt matter. 1:4-6 
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O ,. Difficulties encountered in trying to obtain a closed description of flames are conditioned by the following two crucial 
IJ ' aspects: 

.,— I | (1) Characterized by the flow velocities which are typically well below the speed of sound, deflagration represents 
an essentially nonlocal process, in the sense that the flame-induced gas flows, both up- and downstream, strongly 
affect the flame front structure itself. A seeding role in this reciprocal process is played by the Landau-Darrieus 



(LD) instability of zero thickness flames. 7 ' 8 A very important factor of non-locality of the flame propagation is the 
vorticity production in the flame, which highly complicates the flow structure downstream. In particular, the local 
relation between pressure and velocity fields upstream, expressed by the Bernoulli equation, no longer holds for the 
flow variables downstream. 

(2) Deflagration is a highly nonlinear process which requires an adequate non-perturba- 
tive description of flames with arbitrary values of the flame front slope. As a result of development of the LD-instability, 
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exponentially growing perturbations with arbitrary wavelengths make any initially smooth flame front configuration 
corrugated. Characteristic size of the resulting "cellular" structure is of the order of the cutoff wavelength L c ~ 20Lf 
given by the linear theory of the LD-instability; 3 Lf is the flame front thickness. The exponential growth of unstable 
modes is ultimately suppressed by the nonlinear effects. Since for arbitrary 9 the governing equations do not contain 
small parameters, it is clear that the LD-instability can only be suppressed by the nonlinear effects if the latter are 
not small, and therefore so is the flame front slope. 

The stabilizing role of the nonlinear effects is best illustrated in the case of stationary flame propagation. Numerical 
experiments 9 on 2D flames with = 6 — 8 show that even in very narrow tubes (tube width of the order L c ), typical 
values for the flame front slope are about 1.5 — 2.0. Nonlinearity can be considered small only in the case of small gas 
expansion, 9 — * 1, where one has the O{0 — 1) estimate for the slope, so that it is possible to derive an equation 4-6 
for the flame front position in the framework of the perturbation expansion in powers of (9 — 1). 

This perturbative method gives results in a reasonable agreement with the experiment only for flames with 9 < 3, 
propagating in very narrow tubes (tube width of the order L c ), so that the front slope does not exceed unity. Flames 
of practical importance, however, have 9 up to 10, and propagate in tubes much wider than L c . As a result of 
development of the LD-instability, such flames turn out to be highly curved, which leads to a noticeable increase 
of the flame velocity. In connection with this, a natural question arises whether it is possible to develop a non- 
perturbative approach to the flame dynamics, closed in the sense mentioned above, which would be applicable to 
flames with arbitrary gas expansion. 

A deeper root of this problem is the following dilemma: As was mentioned above, the flame propagation is an 
essentially non-local process; on the other hand, this non-locality itself is determined by the flame front configuration 
and the structure of gas flows near the front, so the question is whether an explicit bulk structure of the flow 
downstream is necessary in deriving an equation for the flame front position. In other words, we look for an approach 
that would provide a closed description of flames more directly, without the need to solve the flow equations explicitly. 

The purpose of this paper is to develop such approach in the stationary case. 

The paper is organized as follows. The flow equations and related results needed in our investigation are displayed 
in Sec. Ill Al A formal integral of the flow equations is derived in Sec. Ill Bl on the basis of the Thomson circulation 
theorem. With the help of this integral, an expression for the vortex mode of the burnt gas flow near the flame front 
is obtained in Sec. IIIII To close the system of equations relating the flow variables on the flame front, potentiality 
of the flow upstream, and of the burnt gas flow after extracting the vortex mode are to be expressed through the 
values of the fuel velocity at the front. This is done in Sec. lIV A II in the form of dispersion relations, using which we 
obtain an equation relating the fuel velocity at the front and the front configuration in Sec. IIV A 21 Complemented 
by the relation defining the local burning rate (the evolution equation), the found equation provides description of 
stationary flames, which is closed in the above-mentioned sense. The developed approach is applied to the particular 
case of zero-thickness flames in Sec. IIVB~T1 Furthermore, it is shown in Sec. IIVB~2l how effects related to the finite 
flame front thickness can be taken into account in the obtained equation. Finally, the case of weak gas expansion 
is considered in Sec. [V] where a single equation for the flame front position is obtained in the third and the fourth 
orders in (9—1). Analytical solutions of the third- and fourth-order equations are found in Sees. IV A~2l and lVB 21 The 
results obtained are discussed in Sec. IVII Appendix lAl contains a consistency check for the expression of the vortex 
mode, derived in Sec. IIIII Some auxiliary mathematical results used in the text are summarized in Appendix iBl 

II. INTEGRAL REPRESENTATION OF FLOW DYNAMICS 

As was mentioned in the point (1) of Introduction, an important factor of the flow non-locality downstream is the 
vorticity production in curved flames, which highly complicates relations between the flow variables. In the presence 
of vorticity, pressure is expressed through the velocity field by an integral relation, its kernel being the Green function 
of the Laplace operator. It should be noted, however, that the jump condition for the pressure across the flame front 
only serves as the boundary condition for determining an appropriate Green function, being useless in other respects. 
Thus, it is convenient to exclude pressure from our consideration from the very beginning. The basis for this is 
provided by the well-known Thomson circulation theorem. Thus, we begin in Sec. Ill Al with the standard formulation 
of the problem of flame propagation, and then construct a formal implicit solution of the flow equations with the help 
of this theorem in Sec. Ill Bl 

A. Flow equations 

Let us consider a 2D stationary flame propagating in an initially uniform premixed ideal fluid. Let the Cartesian 
coordinates (x, y) be chosen so that the y-axis is in the direction of flame propagation, y = — oo being in the fresh 
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fuel. It will be convenient to introduce the following dimensionless variables 

(77, = (x/R, y/R) , (w, u) = (v x /U{ , v y /Uf) , 



IL = (P-P f )/p_U i 
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where U{ is the velocity of a plane flame front, P{ is the initial pressure in the fuel far ahead of the flame front, and R 
is some characteristic length of the problem (e.g., the tube width). The fluid density will be normalized on the fuel 
density /?_. 

As numerical experiments show, stationary flames exist only in sufficiently narrow tubes. Hence, assuming the tube 
walls ideal, and denoting its width by b, we will deal below with spatially 26-periodic flames. More precisely, given 
a flame configuration described by the functions f(rj), w(r), £), u(rj, £), 77 £ [0, +6], where f(j]) denotes the flame front 
position, using the boundary conditions /' = 0, w = for 37 = 0, 6, we continue this solution to the domain rj G [—6, 0] 
according to 

f(v)=f(-v), ™{V, i) = -w(-r],£) , u{t), £) = u(-rj, £) , (I) 

and then periodically continue it to the whole 77-axis. 

In connection with this procedure the following circumstance should be emphasized. All subsequent analysis is 
carried out under assumption that there exists a short wavelength cut-off for the flame perturbations. In other 
words, we consider flames with small but non-zero front thickness. 10 Existence of such a cut-off ensures smoothness 
of the functions under consideration. In particular, it prevents development of singularities such as the edge points 
occurring at the front of a zero-thickness flame, which lead to discontinuities in the values of the flow variables or 
their derivatives. The end points of the flame front, however, still represent a potential source of such discontinuities 
even in the case of flames with non-zero thickness because of the possibility of stream line refraction at these points, 
resulting in a formation of stagnation zones in the flow of burnt matter (see Ref. 11 for more detail). Having imposed 
the boundary condition /'(0) = f'(b) = we thereby exclude this possibility. In view of what has just been said, it 
is natural to assume further that considered as functions of the complex argument, f(rj) together with the "on-shell" 
values of the flow velocity, w{rj 1 /(?))), u(r],f(r])), are analytical functions of 77 in a vicinity of the real axis. This 
assumption, simplifying subsequent analysis, is only technical, and can be weakened if necessary. 

As always, we assume that the process of flame propagation is nearly isobaric. Then the velocity and pressure fields 
obey the following equations in the bulk 



0, (2) 



dvk 1911 1. 1 o ™ 

vi^— = — ^— , k = 1,2, 3 

where (Ci,C2) = (^jOj ( v ii v 2) — (w,u), and summation over repeated indices is implied. 

Acting on Eq. © by the operator Ekid/dQ, where = — em, £12 = +1, and using Eq. one obtains a 2D 
version of the Thomson circulation theorem 

»,| = 0, (4, 

where 

_ du dw 
~ d?] d£ ' 

According to Eq. the vorticity a is conserved along the stream lines. As a simple consequence of this theorem, 
one can find the general solution of the flow equations upstream. Namely, since the flow is potential at £ = —00 
(u = V = const, w — 0, where V is the velocity of the flame in the rest frame of reference of the fuel), it is potential 
for every £ < f(rj). Therefore, 



hoc 

U = 

n — — oc 

w = H{u-u {0) ), (6) 



-f-OO 

J2 u^exp^flnK + iTM?)}, uM = V, (5) 
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where the linear Hilbert operator H is defined by 

H exp{ikrf) = ix(k) exp(ifcr/) , k ^ , (7) 

and 



■ 1 . x>0, 
-1, x < 0. 



In the coordinate representation, H acts on a function a(r]), such that 0(77) — > for 77 — * ±00, according to 



+00 



—00 

It will be shown in the next section how the Thomson theorem can be used to obtain a formal integral of the flow 
equations downstream. 



B. Integration of the flow equations 

Consider the quantity 

, , . _•■ f lnr 

a,i(Q — A Vi = ds v t — , 

s 

where r is the distance from an infinitesimal fluid element ds to the point of observation £, r 2 = — Q) 2 , and 
integration is carried over £ = {r),£ : £ > f(f))}- Taking into account Eq. (|2J, one has for the divergence of a.; : 

-^7- = [ ds v l d l \^— \ = - f ds vSi ( 77— ) = - ( ds di ( v}-^- ) = - [ dh v}-^- , (8) 



dQ J \2tt J J \2tt J J \ 2n J J 2tt 

s s s 

where di = d/dQ, A boundary of E, and dli its element. 
Next, let us calculate £i k d k ei m dia m . Using Eq. ©, we find 

dai f lnr \ d 2 a\ 
dli v t 

\ or] j 

\ A 

d f lnr 
Aai — — / dli Vi 





^ <9a 2 


9ai\ 




— I 


V d?7 







di] J 2n j d£ 2 



drj J 2n 

A 

Analogously, 



d f In r 

£2kdk£imdia m = -Aa 2 - ttz dh Vi—— 

d£, J 2vr 



Together, these two equations can be written as 



In r 

£ikd k £i m dia m = -Aai - di I dl k v k —— 



2n 

Substituting the definition of ai into the latter equation, and integrating by parts gives 



-e ik dk£i m diA 1 v m -d i J dl k v k ^ 

A 

/In t f In r f In r 

dh £i m v m — - di dl k v k — - e lk d k j ds —a . (9) 
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The first two terms on the right of Eq. JSJ) represent the potential component of the fluid velocity, 12 while the third 
corresponds to the vortex component. Our aim below will be to transform the latter to an integral over the flame 
front surface. To this end, we decompose S into elementary ds as follows. 

Let us take a couple of stream lines crossing the flame front at points (77, f(rj)) and (77 + A77, f(rj + A77)) (see Fig.^l. 
Consider the gas elements moving between these lines, which cross the front between the time instants t = and 
t = At. During this time interval, these elements fill a space element So adjacent to the flame front. For sufficiently 
small A77, At, the volume of So 



As 



where 



A77 f'Ar, 
w+At u + At 



dr] 



(u+ - f'w + )Ar]At = v^NA-qAt , 



the subscript "+" means that the corresponding quantity is evaluated just behind the flame front, i.e., for £ = /(fy)+0, 
and ift = Vi+rii is the normal velocity of the burnt gas, m being the unit vector normal to the flame front (pointing 
to the burnt matter). After another time interval of the same duration At, the elements move to a space element Si 
adjacent to So- Since the flow is incompressible, Si is of the same volume as So- Continuing this, the space between 
the two stream lines turns out to be divided into an infinite sequence of S's of the same volume, adjacent to each 
other. Thus, summing over all A77, the third term in Eq. can be written as 



£ -fd k I dl v%v+K{>i.t.fi) 



(10) 



where F denotes the flame front surface (the front line in our 2D case), 

- 00 

K(V,Z,V) = - jjm ^ln{(^-H(,7,nAt)) 2 + (e-S(r),nAt)) 2 } 1/2 Ai 



n=0 



+00 



dtln{( V -R(fj,t)f + ^-E(fj,t)f} 1/2 , 



(11) 



and (H(?7, t), S(ry, t)) trajectory of a particle crossing the point (fj, f(if)) at t = 0. 
Substituting expression illO|) into Eq. gives 

/In 7* f In v f 
dh simVm, — - di I dl k v k-^ T^dk J dl v+cT+K (r], £, 77) . 



(12) 



This representation of the flow velocity downstream will be used in the next section to determine the structure of the 
vortex mode near the flame front. 



III. STRUCTURE OF THE VORTEX MODE 



To determine the flame front dynamics, it is sufficient to know the flow structure near the flame front. As to the 
flow upstream, it is described by Eqs. JSJ, © for all £ < f(rf). Given the solution upstream, velocity components 
of the burnt gas at the flame front can be found from the jump conditions which express the mass and momentum 
conservation across the front. On the other hand, these components are required to be the boundary values (for 
£ = / (77) + 0) of the velocity field satisfying the flow equations. As was shown in the preceding section, the latter can 
be represented in the integral form, Eq. I|12|). Any velocity field can be arbitrarily decomposed into a potential (vf) 
and vortex {vf) modes: 



Our strategy below will be to use the integral representation to determine the near-the-front structure of the vortex 
mode described by the last term in Eq. I|12|) . As to vf, the condition of its potentiality, expressed in the form of a 
"dispersion relation" at the flame front, will eventually close the system of integro-differential equations at the front. 
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Equation i|12|) reveals the following important fact. Up to a potential, the value of the vortex mode at a given 
point (77, £) of the flow downstream is determined by only one point in the range of integration over (fj,t), namely 
that satisfying 



E(fj,t)= V , H(r?,f)=£. 



(13) 



This is, of course, a simple consequence of the Thomson theorem underlying the above derivation of Eq. (|12() . It can 
be verified directly by calculating the curl of the right hand side of Eq. (JTJJl: contracting this equation with ei m d m , 
using EikSim = <5fcm , and taking into account the relation 



Alnr = 2vr(5 (2) (C-C), 



one finds 



+00 



£imd m £ikd k / dl vla + (fj)K(r),£,fj) = dt dl vla + (fj)5(r) - R(fj,t))6(£ - E(fj,t)) 



(14) 



(15) 



Since S(x) = for x 7^ 0, the product of (5-functions picks the point (|13|) out of the whole range of integration in the 
right hand side of Eq. ifHfl) . 

Now, let us take the observation point (£,77) sufficiently close to the flame front, i.e., £ « /(??), [£ > /(»?)]■ In view 
of what has just been said, the vortex component for such points is determined by a contribution coming from the 
integration over 77, t near the flame front, which corresponds to small values of t. Integration over all other 77, i gives 
rise to a potential contribution. 

The small t contribution to the integral kernel -^(77, £,77) can be calculated exactly. For such i's, one can write 



Hfo t) « 77 + w + (fj)t , E% t) « f(fj) + u + (fj)t . 



(16) 



Let the equality of two fields ^1(77,^), ^2(77, U P to a potential field be denoted as tpi = (f2- Then, substituting 
Eq. i|ltj|) into Eq. (|ll|l . and integrating gives 



K(r>,£,fj) = I J *ln{(77-H(77,i)) 2 + (e-S(77,t)) 2 } 


to 

= ~ J dt\n{v\t 2 -2rv + t + r 2 } 1/2 



21 V2 



v+t - 



1 

7T 



dti In \ t{ 



(rv + ) 2 



1/2 



1 



(rv+) 2 



arctan 



(ru|) 2 



-t x + h In <j t\ - (™±L + r 2 



1/2' 



(17) 



Here v + denotes the absolute value of the velocity field at the flame front, and to is assumed small enough to justify 
the approximate equations Ijl6|l . 

As we know, the only point in the range of integration over 77, t that contributes to the vortex mode is the one 
satisfying Eq. ljT3|) or, after integrating over i, 



[77 - 77]tj + (t7) - [£ - f(fj)]w+(fi) = . 



(18) 
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The distance r between this point and the point of observation tends to zero as the latter approaches to the flame 
front surface. Thus, taking [£ — f(rf)} small enough, one can make the ratio t /r as large as desired; therefore, the 
right hand side of Eq. fTTjl is = 



r J I ( rv + 



TTV-l 



2 r 



2 V rv + 



7T 

arcsm 



r^+r^^^+Tic, (i9) 

rv+ rv + 



where "TIC" stands for "Terms Independent of the Coordinates" (??,£)• Denoting 

rv+ 

we finally obtain the following expression for the integral kernel 



K(r].^fi) = — f- + arcsinfi N ) + ftln-} + TIC . 

nv+ I V 2 / e J 



In order to find the vortex mode of the velocity according to Eq. I|12f> we need to calculate derivatives of K. Using 
the relation 



dn _ 1 



- 0- ) , r* = Ci ~ Ci , (20) 



one easily obtains 



M£ i J_ J ( H _ a V J±] ^/2 + arcsin» Vj± 



Equation 121|) can be highly simplified. Consider the quantity 



Ti = - n*±) ""'V 2 + V J± In r . (22) 



Let us evaluate d{Ti . First, we calculate 

T-"— = — " r <3 -^7- = — ' ( 23 ) 



'djTj rj \ dfl v i+ fl 2 
\r v + J \ r 1 r 3 / dri v + 



d ( Vj+_ \ _ v i+ rj _ Q 
v + r 2 r 



(24) 



Second, we note that the vector 



1 r v+ ) JT^W 



satisfies 



PiPi = 1, PiV i+ = , 

i.e., Pi is the unit vector orthogonal to v + . In addition to that, /3j changes its sign at the point defined by Eq. JTHJ. 
Therefore, the derivative of entering T, contains a term with the Dirac (^-function. However, this term is multiplied 
by (arcsinil — 7r/2) which turns into zero together with the argument of the 5-function. Therefore, the product of the 
additional term with (arcsinS! — 7r/2) is to be set zero, in the sense of distributions. Thus, using Eqs. (|20|l . J22J), l|24|) 
one finds 



dTi £1 2 arcsinfi — ir/2 









v + J 



1 



i - n 2 



arcsin — n/2 

(l_f]2)3/2 



on n 



0. 
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We conclude that Tj gives rise to a pure potential. A similar calculation shows that also 



e ik -^^0. (25) 



Therefore, we can rewrite Eq. H21(l as 



dK o 1 ( n n v i+ \ 1 /3., 



0^ = n . (26) 

Finally, substituting this result into Eq. 1)120. noting that the vector £kih is the unit vector parallel to v i+ if 
£iki"iVk+ > 0, and antiparallel in the opposite case, we obtain the following expression for the vortex component 
of the gas velocity downstream near the flame front 

dl X(e pq r p v q+ ) l+ . (27) 

F + 

Having written the exact equality in Eq. I|27|l we take this equation as the definition of the vortex mode. As a useful 
check, it is verified in appendix^ that the obtained expression for v\ satisfies 

(du v /dr)-dw v /dt) + =a- + . 

It remains only to make the following comment in connection with the obtained expression for the vortex mode. 
As is clear from its derivation, Eq. Q27|l is applicable to unbounded as well as bounded flames. In the former case, 
however, the improper integral on the right of this equation is undefined, because integration over the infinite "tails" of 
the flame front around the point satisfying Eq. I|18fl gives rise to a potential contribution which is formally divergent. 
In the case of periodic flames which are of our main concern (see the beginning of Sec. Ill A|) this complication can be 
easily overcome if we specify that the integral is to be understood as 

+2bn 

dr}---= lim f drj--- , (28) 

F -2bn 

so that the contributions of the tails cancel each other exactly. We prefer to work with an infinite ^-interval, rather 
than T] 6 [—b, 6], because of the reasons that will be clear in Sec. llV A"T1 



IV. CLOSED DESCRIPTION OF STATIONARY FLAMES 



After we have determined the near-the-front structure of the vortex component of the gas velocity downstream, we 
can write down a closed system of equations governing the stationary flame propagation. As was explained in the 
Introduction, the term "closed" means that these equations relate only quantities defined on the flame front surface, 
without any reference to the flow dynamics in the bulk. This system consists of the jump conditions for the velocity 
components at the front, and the so-called evolution equation that gives the local normal fuel velocity at the front as 
a function of the front curvature. These equations (except for the evolution equation) are consequences of the mass 
and momentum conservation across the flame front. In Sec. IIV Al we obtain a closed system in the most general 
form, without specifying the form of the jump conditions, and then apply it to the case of zero thickness flames in 
Sec. llV B II 



A. General formulation 



First of all, we need to find the "on-shell" expression for the vortex component, i.e. its limiting form for £ — > /(r;)+0. 
For this purpose we note that in this limit, x( £ ikfiVk+) — * x(v ~ V)i therefore, Eq. 127f) gives 

vl + = jdl X { V -f,) V ^±. (29) 

F 

Let us denote the jump [v i+ (r])—Vi_ (i])] of the gas velocity v j across the flame front as [vi\. Here Vi±(rj) = Vi(rj, f(r])±0). 
Then we can write 
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, r , p , [ ( ~, (g + [v n ])(Vi- + N) ^ 

«*- + [«<] = <+ + y <S %(r? - r))a+ - . (30) 

F 

In every finite order of the asymptotic expansion with respect to the flame front thickness, the jumps [vi] (as well as 
ct+) are quasi- local functionals of the fuel velocity at the flame front, and of the flame front shape (i.e., depend on 
their values and the values of their derivatives of finite order in a given point). Two equations H3()(l . together with 
Eq. © and the evolution equation, v™ = i>"(/), form a system of four equations for the five functions Vi-(r]), vf + (r}), 
and f{rf). To close this system, we need an equation expressing potentiality of the field vf, to be derived in the next 
section. 



1. A dispersion relation for the potential mode 

By the construction, divergence of the last term in Eq. (|12|l . describing the vortex mode, is zero identically in the 
entire downstream region. It is not difficult to see that it retains this property after the number of simplifications 
we have made in the course of derivation of the final expression l|27[) . Note, first of all, that any operating with the 
integral kernel K(rj,t;,fi) itself before differentiation with respect to Q cannot break this property. On the other hand, 
Eq. (|25|1 shows that divergence of the only term, namely Tj, that have been omitted after this differentiation is zero 
identically. In view of the equation divu = 0, one concludes that the potential mode of the velocity field, v p , satisfies 
divv p = downstream, too. Equations rotv p = 0, divv p — allow us to introduce the potential, 0, and the stream 
function, -0, according to 

P _d± _ d± _ dl _ 

U ~ ~ d V ' W - dr, ~ df [ j 

These relations imply that the combination <I> = ip + i<fi is an analytical function of the complex variable z = rj + i£, 
and therefore so is its derivative d^/dz = oj p = u p + iw p . Then, using the Cauchy theorem, we can write 

<JP( Z ) = -L Idz 1 ^- , (32) 
2iri J z — z 
A 

where z e S, and it is assumed that z runs A counterclockwise. In the course of derivation of the expression l|27() for the 
vortex mode, we have been systematically omitting potential contributions. By the construction, these contributions 
are proportional to the integral kernel of the Laplace operator [Cf. Eq. (|12p], and therefore, they generally diverge 
logarithmically at infinity [as, for instance, the last term in the expression (|22[l]. Hence, the integral over the part 
A\_F of the boundary of E in Eq. I|32|l is formally an infinite constant, while the improper integral over F is undefined. 
To avoid appearance of such divergent integrals, we will work below with the velocity derivative duj p /dz, instead of 
u p . Then Eq. is replaced by 

duj p , . 1 f dz dw p . .„„. 

- r {z) = —j>-_ —2. (33) 

dz Zm J z — z dz 

A 

Let us show that this identity can be rewritten as a dispersion relation for dio/dz taken at the flame front. For this 
purpose, let us choose the contour of integration consisting of a large semicircle of radius R, its center being at the 
point zq — + if(0), and of the part of the front F, indented by the circle. Then in the limit z — > F, Eq. I|33|l takes 
the form 



dcjP . . 1 f dz duj p ... 1 7 , dw p . 

— (z+) = —4 — —(z) + - / dtp km — (Re 

dz nij z — z + dz it J R^oo dz 



(34) 



F 

Our choice of the contour implies that the first integral on the right hand side of this equation is understood as 

ZQ-\-2bn 

dz ■ ■ ■ = lim / dz ■ ■ ■ . (35) 

N3rnoo J 
F zo-2bn 
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According to what has been said about behavior of the potential mode at infinity, 



d ^ 

— (Re iv ) ^ as R -> oo for all if G (0, tt). 

Next, differentiating Eq. I|29|l with respect to ?y, taking into account the relation 

dx(x) 



dx 



2S(x) , 



and performing the ^-integration yields 



v w _ Nv^cr+V^ 



(36) 



(37) 



On the other hand, one can write, in view of analyticity of w p , 



V 9 V J + V d£ ) + 



\ dz J 



or, 



\ dz 



(! + </') 



uj = u + iw . 



The right hand side of this equation is explicitly periodic. Hence, the first integral in Eq. Ij34(l is well-defined by the 
rule l|3*5|l. so this equation becomes 



/dujP\ 1 f dz /dujP 



(38) 



A relation similar to Eq. (|38|l can be written for the gas velocity upstream: 



djj 
dz 



1 

7T? 



duj 



Z- \ dz 



(39) 



We have avoided appearance of potentially divergent quantities in the dispersion relation for the gas flow downstream 
at the cost of increasing its differential order by one. 13 Otherwise, we would have had to carry out all intermediate 
calculations for a finite domain E, and to require all divergences to cancel in the final equation for the flame front 
position in the limit of infinite front length. In both cases, the finite constant in this equation (an integration constant 
in the first case, or a finite remainder after the cancellation in the second) is fixed by the normalization condition 



V = 



+b 



drj N. 



(40) 



2. The integro-differential relation between u>- and f. 

We are now in a position to write down the main integro-differential equation relating the values of the gas velocities 
at the flame front with the flame front position. To this end, we first differentiate Eq. (130(1 with respect to 77 doing 
the integral as before, and rewrite the result in the complex form: 

(w_) +H =(w£) +Na + ^ - ^ 2]) . (41) 



Acting on this equation by the operator 



mj z — z+ 
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taking into account Eqs. Q38[l . I|39[) . and the relations 



gives 

2 <„_)< + (i + rfc) {h- - ^ ( ° ;+ ( ^ ) +M) } - o , <m 

where 

£ = l + if f dfj 



IT J Z — Z + 
F 

Written longhand, the action of the operator 3~C on a function a(r]) is 



+ 0O 

n w = i±m 7 - • (43) 



J{ has properties similar to the Hilbert operator H. In particular, it is shown in appendix IB II that 

IK 2 = -1. (44) 
In terms of !K, cj_, the dispersion relation (|39|l takes a more compact form 

(l - ift) = . (45) 

In view of the property (|44|) . this relation is fulfilled by any w_ satisfying Eq. I|42|). 

Equation l|42|l is the main integro-differential relation between the fuel velocity at the flame front, and the flame 
front position. It can also be rewritten in terms of lu + as 



= . (46) 



Together with the jump conditions and the evolution equation which has the general form 

i£ = l + S(u_, w-,f), (47) 

where S is a quasi-local functional of its arguments, proportional to the flame front thickness, Eq. I|42|) provides a 
closed description of stationary flames. 

B. Equation (1421) in lowest orders of the e-expansion 

In this section, the general results obtained in the preceding section will be applied to thin flames in the first two 
orders of the asymptotic expansion with respect to the flame front thickness e = L{ / R. 

1. Zero-thickness flames 
For zero-thickness flames, the jump conditions for the velocity components have the form 

M = ^ , (48) 
M = -f' 9 -^' (49) 



12 



and hence, 

H = {e-i) l —^f- , (50) 

while the evolution equation 

v n _ = 1 . (51) 
We see that the jumps are velocity- independent, and 5 = 0. Also, it follows from these equations that 

[v n ] = 9-1, [v 2 ] = e 2 -i. 

It remains only to calculate the value of the vorticity at the flame front, as a function of the fuel velocity. This can 
be done 2 directly using the flow equations With the help of Eqs. (5.32) and (6.15) of Ref. 2, the jump of the 

vorticity across the flame front can be written, in the 2D stationary case, as 



° - ( «- (53) 



where 

f) = I in 4- 

N J rfyy 

Differentiating the evolution equation written in the form 

u_ - f'w- = N , (54) 
and writing Eq. I|52|) longhand, expression in the brackets can be considerably simplified 

Dw- + f Du_ + —Df = w_w- + — — h — — h / u_w_ + — 

(w 2 _)' (u--N)' (N 2 -l)u'_ N' 
= + N + ~ j^— + «-(«- - N) + — = «'_«_ + w'_w^ . (55) 

Since the flow is potential upstream, we obtain the following expression for the vorticity just behind the flame front 14 

a+ = — w~^ u '~ u ~ + w '- w -} ■ ^ 

Substituting these expressions into Eq. (|4*2"|) yields 

("-)' + ^ (l + ^) | + \ («- + (* " V^f) [ ln d--| 2 + 92 !)]'} = • 

(57) 

Equations 154|l and (|57l) provide the closed description of stationary zero-thickness flames in the most convenient 
form. Account of the heat conduction - species diffusion processes inside the thin flame front changes the right hand 
side of this equation to 0(e). This modification will be considered in the next section. 



2. Account of the transport processes in the linear approximation 

An equation relating lo_ with / can be obtained also in the case of flames of nonzero thickness following the lines 
of the above derivation of Eq. JS7J. For instance, in the first order in the flame front thickness e, one has to use in the 
general Eq. (|42|) the jump conditions and the evolution equation derived in Ref. 2. The resulting equation turns out 
to be much more complicated than Eq. (|57|l . However, the main purpose of taking into account the inner structure 
of the flame front is to provide a short wavelength cutoff for unstable flame perturbations, which is necessary for the 
very existence of stationary configurations of curved flames. On the other hand, for many purposes it is sufficient 
to consider the transport processes in the linear approximation, while full account of the nonlinear coupling of these 
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processes to flame hydrodynamics is of minor importance in practice. Therefore, we can generalize Eq. (I57|l to the 
case of flames of nonzero thickness almost without calculation using the results of the linear theory of flame front 
instability. In the linear approximation, cj_ = 1, hence, the left hand side of Eq. 1)57(1 reduces to 



H-i)f 



8-1 
2 

On the other hand, omitting the time derivatives in the linear equation 3 describing evolution of the front perturbations 
gives formally 

Hf> = ^f" 

where A c = L c /R = 0(e) is the dimensionless cutoff wavelength. We conclude that the desired generalization of 
Eq. JS7JI reads 15 

+ «zl (i + dt) { + i + (o - wr-p +«■-!)]'} 

9-lXr. 



i + m)f". (58) 



2 2tt 

Equivalently, the main system of equations can be rewritten in terms of W+ , u + : 

^ + « - + ' + - (1 + - £ r ) _ o , 

u+ - /'io+ = 6W . 

We have written Oi instead of H in the e-terms in order to formally preserve the dispersion relation 1(45(1 . However, 
it should be kept in mind that by the construction, the e-term is to be treated linearized. 

Most probably, Eq. 1(58(1 can be solved in full generality only numerically. However, it admits analytical solutions 
in lowest orders of the asymptotic expansion for — > 1, considered below. 



V. THE SMALL (9 - 1) EXPANSION 

Although our main concern in this paper is the non-perturbative description of flames with arbitrary 8, it is of some 
interest to apply the above results to the case of weak gas expansion. One reason is that this case allows considerable 
simplification of Eqs. 1(54(1 . ((58(1 which can be reduced to a single equation for the flame position. Another is that 
this is a good place to illustrate the role of the relation l(4l)J) in our approach. The well-known subtle point of the 
analytical theory of nonlinear flame propagation is to ensure that the constant V and the function f(j]), which play 
the role of the eigenvalue and the eigenfunction of the nonlinear equation for the front position, respectively, satisfy 
this obvious condition. 

In the present formulation, this problem does not arise at all. To see this, it is sufficient to note that the constant 
V does not appear at all either in Eq. ((58(1 . or in the evolution equation 1(54(1 . Therefore, Eq. ((40(1 is to be considered 
simply as a relation that allows one to express the flame velocity through the constant of integration C [Cf. discussion 
in the end of Sec. IIV A 1( . In lowest orders of the small (6 — 1) expansion, this connection between V and C can be 
followed out in detail. It will be shown in Sec. IV Al that in the first post-Sivashinsky approximation, Eqs. ((54(1 . I(58|) 
reduce to the Cambray-Joulin 16 version of the well-known Sivashinsky-Clavin equation. 4 In Sec. lVBl an equation for 
the flame front position will be obtained in the second post-Sivashinsky approximation, which represents a corrected 
version of the equation obtained by Kazakov and Liberman 5-6 (called there the fourth order equation). 



A. The first post-Sivashinsky approximation 

1. Equation for the flame front position 



The form of the £-term in Eq. 1(58(1 shows that the cutoff wavelength A c ~ 1/a, a = 8 — 1. This implies that 
/' = 0(a), and that differentiation of a flow variable with respect to rj increases its order by one. Furthermore, it 
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follows from Eq. l{5~4*|) that u_ = 1 + (3(a 2 ), since V = 0(a 2 ). Then Eq. © tells us that also W- = 0(a 2 ). To carry 
out expansion of Eq. I|58|) . we need also to find the asymptotic action of the operator IH. A general procedure of 
consistent asymptotic treatment of this operator is formulated in appendix IB 21 At the third order in a, according to 
Eq. (|ljl2|) . OH acts on the expression in the curly brackets on the left hand side of Eq. I|58|) as the Hilbert operator. 
Thus, within accuracy of the fourth order, Eq. (|58|l can be rewritten as 

( "- } + — \ l + lH ) {-^ + ^9--2^ f ) =°- 

Expanding the last term in this equation, noting that with the required accuracy |w~| 2 in the parentheses can be 
taken to be 

M 2 = i + (/') 2 , (60) 

and integrating gives 

w _-^i(i + ^)( i /' + ^r) = c ) (6i) 

where C = C\ +1C2 is the constant of integration. Extracting the real and imaginary parts of this equation, we obtain 

"- - V(-* / ' + s / ")= c " (62) 

■»- " V('' + 5^")= c " <63 > 

Integrating Eq. (|63|l over [—6, +b], and taking into account periodicity of the function f{rj) shows that C2 = 0, because 

+b 

J dri W-(rj) = 

-b 

in view of the definition 0). Next, multiplying Eq. Q63JI by /', subtracting it from Eq. I|62|l . and using the evolution 
equation l|54(l . we arrive at the single equation for the flame front position 17 

f (/f = ^ (-*/•+ &)+a-i. (64) 

The constant C\ can be expressed through the flame velocity V using the normalization condition i|40|) . Namely, 
integrating Eq. I|64|l over [—6, +6], and applying this condition gives 

9 (V - 1) = Ci - 1 , 

thus bringing the third order equation to the form 

-9W+ 6 - {f<f = ^ (-#/' + ^/") , (65) 

where = V — 1 is the flame velocity increase due to the front curvature. This equation exactly coincides with 
the stationary version of the corrected Sivashinsky-Clavin equation (obtained by setting ^ = — W, and restoring the 
e-term in Eq. (10) of Ref. 16). 



2. Solution of the equation \65\l 

The third order equation (|65J) is of the same functional structure as the Sivashinsky equation. 1 Therefore, it can be 
solved analytically using the method of pole decomposition. 18,19 We look for a 2&-periodic solution of Eq. (|65|l in the 
form 



2P 



k=l 



(66) 
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The amplitude A and the complex poles T)k, k = 1, ...,2P are to be determined substituting this anzats into Eq. I|65[l . 
Since the function f(rj) is real for real rj, the poles come in conjugate pairs; P is the number of the pole pairs. 
Using the formulas (see Refs. 18,19 for more detail) 



fc=l 

cot x cot y = — 1 + cot(ir — y)(cot y — cot x) , 
it is easily verified that Eq. is satisfied by f(rj) taken in the form of Eq. provided that 

A c 0-1 



2P 



1 + i sign(Im r] k ) cot — (rj 



>]} 



sign(x) = — , 

a; 



A = — 



W 



2tt (9 



l) 2 PXr 



1 



PA, 



and the poles k = 1, 



20 2 26 

, 2P , satisfy the following set of equations 
\ 2P 

% L COt oiifa 



26 



i sign(Im 7^) 



rc = l 



0, fe = 1, 



,2P. 



(67) 



(68) 



It is seen from Eq. 1(67(1 that the found solution Ij66(l is not unique: there is a number of solutions corresponding to 
different numbers P of poles. To identify the physical ones, a stability analysis of the solutions is required which, 
of course, cannot be carried out in the framework of the stationary theory. However, as we have mentioned above, 
the functional structure of Eq. (|65|l is very similar to that of the stationary Sivashinsky equation. Under assumption 
that the non-stationary versions of these equations are also similar, the stability analysis performed in Ref. 20 will 
be carried over the present case. According to this analysis, for a given (sufficiently small) b, there is only one stable 
solution. This solution is identified as that maximizing the flame velocity. In addition to that, the poles of the stable 
solution are required to be vertically aligned in the complex 77-plane. For such a "coalescent" solution, a simple upper 
bound on the number of poles can be obtained from Eq. 168(1 . Namely, for k — fco with r)k uppermost, one has 



1 = ^ , 
2b ^ 



2P 

^ coth 



^(Im^o-Im^)] > ^(2P-1). 



(The equality holds, if Im r]k = 00) Then Eq. I|67|l tells us that the maximum W is reached for 



P, 



Int 



Ac 



Int(x) denoting the integer part of x. Thus, the flame velocity increase W s of the stable solution takes the form 

PraaxA/ 



2h 



1 



2b 



(69) 



where 



(70) 



B. The second post-Sivashinsky approximation 

Before going into details of derivation of the correct fourth-order equation, let us take a pause to identify the origin 
of the failure of equations derived in Ref. 6 to satisfy the condition (|40|l . It is traced to the choice of the constant 
term in the decomposition Eq. (19)' of the potential mode 21 of the downstream flow (the formulas cited from Ref. 6 
are distinguished by the prime). In Ref. 6, this term is taken to be equal to QV. It reappears later in the right hand 
side of Eq. (39)' which is an integral of the main Eq. (40)'. It was proved 6 that Eq. (40)' is valid up to terms of 
the order a 6 . However, this does not mean that Eq. (39)' is valid up to terms of the order a 5 . The point is that the 
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value 8V of the constant term in the potential mode is accurate only up to 0(a 3 )-terms. In other words, the proper 
decomposition of u p should read 

+00 

Up = 8V + P + u p , u p = I dkuk exp(— |fc|£ + ikrj) , 



where = 0(a 3 ) is some constant, instead of Eq. (19)'. Correspondingly, Eq. (39)' is to be substituted by the following 

6 = -H\n- 



8V + /3 

Since 9, V, fi are all O(l), this implies that an additive constant of the order 0(a 3 ) is missing in Eq. (39)'. The value 
of the constant (3 is fixed eventually by the condition (|40|l . As in the preceding section, it can be expressed through 
W explicitly by averaging the resulting equation for the flame front position. It is not difficult to verify that following 
this way, one obtains exactly Eq. ll6*B|) instead of Eq. (70)'. 



1. Equation for the flame front position 



Turning to the derivation of the fourth order equation, we write 



|w _ |9 = l + U'n^+w-f = j + (f)2 + 2fw _ + Q(a4) ( (n) 



instead of Eq. i|60J) . Substituting this into Eq. 1)58(1. using Eq. (jB 10|> - and integrating yields 

N ' 26> 2n~ 




f"(fj)im-m]=c, (72) 

which replaces Eq. I|61l) . Next, substituting the third order result for u>_ in the curly brackets in this equation, one 
finds with the required accuracy 

+00 

- - V + -) + - '-wrf 5f;«* m '>> - - c ■ 

—00 

When extracting the real and imaginary parts, one notes that the last term in the latter equation is purely imaginary, 
and hence it affects only the W- component of the fuel velocity. Since this term is of the fourth order, it can be 
omitted because W- is multiplied by /' in the evolution equation. Thus, we have 



2 \N J 2 2tt 

while for w_ one can still use the third order Eq. I|63(l (with C2 =0). Apart from the e-terms which are treated here 
differently from Ref. 6, equation (|73|l is nothing but the expression Eq. (64)' for v (the ^-independent counterpart 
of it, introduced in Ref. 6), in which the constant term is modified as discussed above. Substituting from 
Eqs. (|63|l . (|73J) into the evolution equation Q54JI gives 

if ( il , tt\ ( J \ , hi 



Averaging this equation along the flame front as before, noting that f'(—rj) = —/'(»?), as a consequence of Eq. (Q, 
and applying the condition Eq. 14U|) one can express the constant C\ through W. Thus, we obtain 

+ (^) 

As in the case of the third order equations (|65|l and (70)', the fourth order equations l|74l) and (73)' coincide exactly 
upon differentiation. To see this, one has to expand the iV-factors in Eq. (|74|l . to linearize the e-terms in Eq. (73)', and 
to take into account that in the course of derivation of this equation, the Sivashinsky equation was used to transform 
some of the nonlinear terms. 
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2. Solution of the equation [74t 



Remarkably, Eq. I|74[l can be further simplified, and brought to the form similar to that of Eq. I|65|) . This can 
be done by expressing the terms of the third and fourth orders in /', which appear upon expanding the N- factors, 
through the lower-order terms using the Sivashinsky equation. Namely, rewriting this equation as 



(/') 



W 



1 



Hf + 0(e) 



2 2 

and taking its square, the term (f) 4 /8 coming from the N on the left of Eq. (|74(1 is substituted by 



(75) 



W 2 — W{9 — l)Hf 



which is already of the second order in /'. The e-terms have been omitted because we neglect all nonlinearities related 
to the transport processes. Next, multiplying Eq. H75fl by /', and acting by the Hilbert operator, the term H(f') 3 /2 
appearing on the right of Eq. (|74|l upon expanding H(f/N) becomes 



H 



(/ 



r\3 



WHf 



- 1 



Hf) -(/')' 



where we have used the well-known identity 



2H (aHa^) = (hoX 



After substitution of these expressions, Eq. (|74() can be rewritten in the following form, within accuracy of the fourth 
order, 



20 2 (ff 
(0 + l) 2 



ew- 



w 2 



i 



-Hf + ^f" 

ITT 



(76) 



Proceeding as in Sec. IV A~2l we look for a periodic solution of Eq. (|76[1 in the form Eq. (|66[) . and find 

A, 0-1 



ew 



A = — 
W 2 



2ir AO 2 y 

> 2 - I) 2 PX r 



1) 



2b 



1 - 



2b 



while equations for the pole positions remain the same [see Eqs. (|68|) ]. As before, we identify the stable solution as 
that maximizing the flame velocity. In particular, the maximum of the flame velocity increase for a given 6 is 



VKmax = -t 



160 2 



Because of the factor 1/16, the second term under the square root represents a relatively small correction to the first 
even for values of not close to unity, so the above expression can be simplified to 

^ max = L__L. (77) 

Figure [3 compares the dependencies of the maximal flame velocity increase on the gas expansion coefficient, given 
by the Sivashinsky equation, and its corrections, Eqs. I|65|l . I|76(l . with the experimental data. 9 ' 22 It is seen that the 
theoretical curve W max (0) approaches the experimental marks as we pass from the Sivashinsky equation to the more 
accurate Eqs. I|65|) . I|76|l . However, one should remember that this improvement can be trusted only for sufficiently 
small values of (9— 1). For large values of the gas expansion coefficient it is rather a lucky accident. Only integration 
of the exact Eq. can give reliable results for the practically important values of 9 . 
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VI. DISCUSSION AND CONCLUSIONS 



The results of Sec. lIVI solve the problem of closed description of stationary flames. The complex Eq. (|42|l determines 
the on-shell distribution of the fuel velocity as a functional of the flame front configuration in the most general form, 
while the evolution equation plays the role of a consistency condition which gives an equation for the front position 
itself. We have shown, furthermore, that in the case of zero-thickness flames, the main Eq. (|42() takes the form JS7J. 
This equation is universal in that any surface of discontinuity, propagating in an ideal incompressible fluid, is described 
by Eq. (|57|l whatever internal structure of this discontinuity be. The latter shows itself in the 0(£)-corrections to 
this equation, where e is the relative thickness of the discontinuity. A simple comparison with the results of the 
linear theory has shown that the linear account of the heat conduction - species diffusion processes in the flame front 
modifies Eq. JHTJ to Eq. fifty. 

Next, some technical remarks are in order. It is difficult to say at the moment whether Eq. (|58|l admits further 
simplification. What can be said, on the other hand, is that its analytical structure does not present serious problems 
for numerical integration. Indeed, the keystone of this structure is the integral operator !K. But its properties are 
much like those of the Hilbert operator H which is well-known how to deal with. From the theoretical point of view, 
Eq. H58|) is convenient for constructing various approximate descriptions of stationary flames. In particular, it is well 
suited for developing the small {9 — 1) expansion which we have carried out in Sec. [V] Specifically, it was verified 
in Sec. IV Bl that at the second post-Sivashinsky approximation, the exact equations derived in the present paper 
reproduce the fourth order equation obtained in Ref. 6, up to an additive constant. This agrees completely with the 
main result of Sec. IIIB of Ref. 6, according to which Eq. (40)' correctly approximates the exact equation for the 
flame front position up to terms of the sixth order in (8 — 1). The difference in the values of the additive constant 
was found in Sec. IV Bl to be the result of an improper integration of Eq. (40)', namely of an incorrect separation of a 
constant term in the Fourier decomposition of the potential mode of the flow velocity downstream. 

The results presented in this paper resolve the dilemma stated in the Introduction in the case of 2D stationary flames. 
Since in our considerations we have extensively used specifically 2D mathematical tools, the question of principle is 
whether the results of this paper can be carried over to the 3D case, and further to the general non-stationary case. 

APPENDIX A: CONSISTENCY CHECK FOR EQ. (1271) 

After a lengthy calculation in Sec. IIIII we obtained the following simple expression for the vorticity mode near the 
flame front 

^ v r ! r a + Vi + 

4 



< = I dl x(e Pq r p v q+ ) - . (Al) 



As this important formula plays the central role in our investigation, a simple consistency check will be performed 
here, namely, we will verify that vf given by Eq. ijAlfl satisfies 

[du v /dr] - dw v /di) + = a+ . (A2) 

Contracting Eq. IjAlfl with Skidk, and using relation (|36[> . one finds 

£kid k Vi = dl 8(£ pq r p v q+ )e ki £ km v m+ +a+ 2 ' + = J dl 6(e pq r. p v q+ )v™(r + . (A3) 



The argument of the (5-function turns into zero when the vectors and Vi+ are parallel. Near this point, one can 
write 

e pq r p v q+ w rv + (j) , 

where </> is the angle between the two vectors. On the other hand, the line element, dl, near the same point can be 



written as 



dl = — — -defy — rd§— 
sin w v. 



as a simple geometric consideration shows, see Fig. [3] 
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Substituting these expressions into Eq. i|A3|) . and taking into account relation 



S(ax) = -. — tS(x) 
\a\ 



we finally arrive at the desired identity 



v 



EkidkVi = I # *(">+0)«+o+ = / # = CT + • (A4) 



+ 

F F 



It should be noted in this respect that the identity Eq. i|A2(l is only a necessary condition imposed on the field w". 
Playing the role of a "boundary condition" for the vortex mode, Eq. (|A2I) is satisfied by infinitely many essentially 
different fields, i.e., fields which are not equal up to a potential. It is not difficult to verify, for instance, that the 
velocity field defined by 

Vi = dl x{£pqrpn q )^-^- 



also satisfies Eq. (|A2|I . and the difference v\ — is essentially non-zero. 

By the construction of Sec. II I II given by Eq. (|A1(I is essentially the only field that satisfies the flow equations 
©-©■ 

APPENDIX B: SOME PROPERTIES OF THE OPERATOR % 
1. Proof of the identity (I44B 

In the course of derivation of Eq. I|42(l , we have introduced an operator IK defined on functionals of the flow variables 

by 

+ 00 

v j 7r j v-v + ntiv) - fwl 

—00 

and used repeatedly the following important identity it satisfies 

JC 2 = -1 . 

To prove this identity, it is convenient to represent the right hand side of Eq. (|B1I) as an integral over the contour 
C\ = Cf U in the complex 77-plane, shown in Fig. ^ 

* o)w = i±i£M/ d , ; |gI, (B2) 

Ci 

where z = rj + if(rj), z = fj + if(fj), and C\ is chosen such that all singularities of the integrand, except the pole at 
fj = r], remain above Cy, or below C-j". Under our assumption about analytical properties of the functions involved 
[see discussion below Eq. Q], such contour always exists. Now !K 2 takes the form 

! + */'(»?) f f Jm afa) 

z\ — z 



(&a) (,) = ^™ / ^ / dm , ,1 = Vl + i/fo) , (B3) 



C\ C2 



where the contour C2 = C 2 U of integration over 771 comprises C\ (see Fig.0J. Changing the order of integration 
in Eq. (|B3Jl . using the formula 
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and taking into account that the logarithm gives rise to a nonzero contribution only if the arguments of the functions 
z — z and z — Z\ run in opposite directions when z runs the contours , we obtain 



l + if{rj) f dz 



4tt 2 



drji 



a( m ) l + if(v) 



z\ — z 



2iri / drji 



C-2 



[1 + if'(v)}^- I dm = -[1 + if '{ V )] ■ res 

Z7T J Z\ — Z 



C 




Z\ — z 



(vi-v) + i[f(vi)-fin)} 



-a{r)) 



and thus 



(lK 2 a) (t?) = -a{r,) 



which was to be proved. 



2. Asymptotic form of IK for 9 — > 1 

It was shown in Sec.|V]that the exact equations derived in Sec. HVAl constitute a general framework for perturbative 
treatment of flames within the small (9 — l)-expansion. An important step of this program is the construction of an 
asymptotic expansion of the operator !K. 

Let us consider first an IH-transform of a total derivative (this is the case we dealt with in Sec. Iv|) 

M ^ - 1 - ±im -f dfj ~ + a !«\ • (B4) 

V J 7T J Tj - r\ + i[f{ri) - f{r))\ 

— oo 

The aim of the subsequent transformations will be to develop an expansion of this integral in powers of [f(fj) — f(rj)] 
which is small for all fj, rather than in powers of (77-/7). For this purpose we first rewrite it as an integral over the 
contour C x [Cf. Eq. (1521 1 



Integrating by parts then yields 



(&»') (77) 



(&»') (ji) 



if in) f Mv) 



2tt 

Ci 



j_d f d ~^mi t (B6) 

2tt drj J z — z 



where C\ denotes a contour in the complex z-plane, which is run by z = 77 + 1/(77) when 77 runs C\. Next, let us make 
the following change of the integration variable in Eq. 1B6I) 



z -> z = z - 7/(77) • 

This gives 



J_d f^afm (B7) 

2tt drj J z — 77 



where the contour C\ of z-integration is obtained by a uniform vertical shift of C± by an amount f(rj). Note that 
z = 77 + i[f(fj) — fiv)] — V + 0(a). Hence, within the framework of the asymptotic expansion in a, C\ is to be 
considered as a small deformation of the initial contour C\. Taking a small enough we can always secure C\ from 
crossing singularities of the integrand (including the pole 77 = 77) under this deformation. Assuming this, we deform 
C± back to the real axis, and obtain 

+00 

(^a')(r ] )= 1 -±-f dfj^m, (B8) 
V / tt drj J rj — rj 
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where 77(77) is a function of the real variable 77, defined by 

V = V + i[f(fj) - f{rf)}. 
Resolving this relation with respect to 77 in the form of a series 

V = v-i[f(f}) - f(v)} 

= v - mv) - m] rmiiv) - fm + \fmfif,) - im 2 +■■■ 
= n- i[m - m] ~ rmm ~ m] + * wm 2 - im 
+lf"mm)-m?+-- , m 

and substituting it into Eq. (|_B8fl . we finally arrive at the following expansion of IKa' 

(fta') (77) = ~ j A {a(77) - 7a'(77)[/(77) - f(r>)} - a' (fj) f (fj)[f(fj) - f(r,)} 

— OO 

+ia'(v) [f'm 2 [f(fj) - f(v)} \a"{ft)\m - fin)? + •••}• (BIO) 

Since each term in the curly brackets in the integrand on the right of Eq. (fir? 1 0|) is periodic, the corresponding improper 
integrals are well-defined by the rule Ij35(l . 

Now it is straightforward to write down the result of asymptotic action of J{ on a general 0(77). Namely, using 
d/di]H = Hd/drj in the first term in Eq. (|B10j) . and making the substitution 0/(77) — > a(rj) throughout this equation 
yields 

+00 

(%a) (77) = (Ha) (77) + ~ j ^ {-ia(fj)[f(fj) - f(t,)] - a(fj)f'mm) ~ M] 

— OO 

+ia(fj) [f'm 2 [f(fj) - f{ v )] - \a'{f))[f{n) ~ f(v)? + •••}• (BH) 

In the course of derivation of the third order equation in Sec.[V] we found it necessary to determine the asymptotic 
action of IK on the curly bracket in the left hand side of Eq. (|58|l taking into account 0(a 3 ) terms. In this case, 

a' = -if" + u'_=0(a 2 ), 

and therefore, it is sufficient to keep only the first term in the integrand of Eq. HljlOfl . which gives immediately 

(Ad') (77) = (Ho!^ (77) + 0(a 4 ) . (B12) 
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Figure captions 
Fig.l: Elementary decomposition of the flow downstream. 

Fig. 2: Dependence of the maximal flame velocity increase on the gas expansion coefficient, given by the Sivashinsky 
equation (dotted line), and its corrections - Eq. (|65fl (dashed line), and Eq. (|76[) (full line). The marks represent the 
results of numerical experiment 9 ' 22 (accuracy of the experimental data is about f0%). 

Fig. 3: Near-the-front structure of the flow downstream, = nv\ is the normal component of the velocity. Since 
the observation point (?7,£) is close to the flame front, the stream line and the part of the front near this point can 
be considered straight. 

Fig. 4: Contours of integration in Eqs. (|B2|I - (|B7|) 
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